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I. INTRODUCTION 



It becomes urgent to study the generation and the propagation of gravitational waves 
in detail, because of the increasing expectation of direct detection of gravitational waves by 
Kilometer-size interferometric gravitational wave detectors, such as LIGO [|I|] and VIRGO 
0] now under construction. Coalescing binary neutron stars are the most promising can- 
didate of sources of gravitational waves for such detectors. Gravitational waves generated 
by them bring us informations not only on various physical parameters of neutron star 
||.[|] , but also on the cosmological parameters [^-fl if and only if we can make a detailed 
comparison between the observed signal with theoretical prediction during the epoch of 
the so-called inspiraling phase where the orbital separation is much larger than the radius 
of component stars |J. The problem is that, in order to make any successful comparison 
between theory and observation, we need to know the detailed waveform generated by 
the motion up to 4PN order [Hj which is of order e 8 higher than the Newtonian order, 
where e = orbital velocity/light speed. 

Blanchet and Damour have developed a systematic scheme to calculate the waveform 
at higher orders, where the post-Minkowskian approximation is used to construct the ex- 
ternal field and the post-Newtonian approximation is used to construct the field near the 
material source [jTT| |3~3]| . In the post-Minkowskian approximation, the background geome- 
try is the Minkowski spacetime where linearized gravitational waves propagate |l^JTT|JT2]| . 
The corrections to propagation of gravitational waves can be taken into account by per- 
forming the post-Minkowskian approximation up to higher orders. In fact, Blanchet and 
Damour obtained the tail term of gravitational waves as the integral over the past history 



of the source p2| , pT3| . They introduced a complex parameter B ||11|| , and used the analytic 
continuation as a useful mathematical device to evaluate the so-called log term in the tail 



contribution |[12|JT3[| . The method is powerful, but it is not easy to see the origin of the 
tail term. Will and Wiseman [|15[] have also obtained the tail term associated with the 



mass quadrupole moment, by generalizing the Epstein- Wagoner formalism [IIS]. In the 
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following, we shall study the waveform formula from different point of view because of 
its importance. Namely we shall use the Green's function for the wave operator in the 
Coulomb type potential generated by the mass of the source. This method is straight- 
forward and enables us to see explicitly how the tail term originates from the difference 
between the flat light cone and the true one, which is due to the mass of the source 



GMj c 2 as the lowest order correction. Schafer JlTj and Nakamura [18| have obtained the 
quadrupole energy loss formula including the contribution from the tail term by studying 
the wave propagation in the Coulomb type potential. 

This paper is organized as follows. In section 2, we derive the formula for the tail 
of gravitational waves by using the Green's function in the Coulomb type potential. In 
section 3, we discuss its relation with previous works. We use the units of c = G = 1. 

II. GRAVITATIONAL WAVES IN COULOMB TYPE POTENTIAL 

We wish to clarify that the main part of the tail term is due to the propagation of 
gravitational waves on the light cone which deviates slightly from the flat light cone owing 
to the mass of the source. For this purpose, we shall work in the harmonic coordinate since 
the deviation may be easily seen in the reduced Einstein's equation in this coordinate. 

_ h a/3 )h^ af3 = -lGKQf" + h%h^ a , (2.1) 

where 7]^ is the Minkowskian metric, h fiu = rj^ — \/—gg^ u and 

e^ = H7)(T^ + C L ). (2.2) 



Here, t^ L is the Landau-Lifshitz pseudotensor Since we look for 1/r part of the 

solution, the spatial derivative is not relevant in the differential operator. Therefore, it is 
more convenient for our purpose to transform Eq.( [2.1|) into the following form 

(□ - h 00 d do)h^ = -16nS^, (2.3) 

where □ = rj^ u d^d v and we defined as 



s »v = qv» _ J_ {h^ J^P + 2h 0i h^ oi + li ; 'ir ;i ). (2.4) 
Substituting the lowest order expression for h 00 , we obtain our basic equation 

U M hr = f»\ (2.5) 
where the effective source f^ u is defined as 

=S^--^-( h 00 - —) h^, (2.6) 



16tt 

Here, M is the total mass of the source and O m is defined as 

4Mx d 2 



r ' df z 

where A is the Laplacian in the flat space. In case of the source with strong internal 
gravity like binary neutron stars, M may be regarded as the ADM mass pl |, since 1/r 



part of h 00 in the harmonic coordinate becomes the ADM mass. In this way, our derivation 
may be applied to such compact sources. 

The solution for fl2.5|) may be written down by using the retarded Green's function in 
the following form 

hT{x) = J d 4 x'G^\x,x')f^(x r ). (2.8) 

The retarded Green's function is defined as satisfying the equation; 

n M G { +\x,y)=5 A (x-y), (2.9) 

with an appropriate boundary condition. 

The Green's function satisfying Eq.fl2.9|) can be constructed by using the homogeneous 
solutions for the equation 

□m# = 0. (2.10) 
The homogeneous solution for Eq. fl2.10D takes a form of 



e-^fiWYUO^), (2.11) 

where we defined 

p = ur. (2.12) 
Then the radial function fi{p) = pfi(p) satisfies 

so that Eq.( |2.11| ) is a solution for Eq.( |2.10|) . Thus we can obtain homogeneous solutions for 
Eq.( |2.10| ) by choosing fi(p) as one of spherical Coulomb functions; uf^{p\ 7) and Fi(p; 7) 
with 7 = —2Mu. Here, we adopted the following definition of the spherical Coulomb 
function PH as 



F{ P] 7) = c ie ip p l+1 F(l + 1 + z 7 |2Z + 2| - 2ip), 
u, (±) (p; 7) = ±2ie Tiai c l e ±ip p l+l W 1 (l + 1 ± ij\2l + 2\t 2ip), (2.14) 



where q and <j; are defined as 



ci 



|r(/ + i + i 7 )| 



(2/+1)! ' 

<Ti = ar g r(/ + l + z 7 ). (2.15) 

Here, F and Wi are the confluent hypergeometric function and the Whittaker's function 
respectively. These spherical Coulomb functions have asymptotic behavior as 



(±) 

u\ ~ exp 



±z(p — 7 ln2p — f° r r ~~ * 00 ' (2.16) 
and 

F/ ~ Qp' +1 for r -> 0. (2.17) 
Thus we obtain the retarded Green's function in the following form. 
G ( m\x,x') = $> iCTi J dusgn{iv)(y +ujlm {x)y Sujlm *{x')6{r - r') 



tin 



+ ty Sujlm (x)ty +ujlm *(x')6(r' - rj), (2.18) 



where we defined ^f +ulm (x) and \^ Swlm ( x ) as 



+oj Imi 




{x) = \^e-^p- 1 u^\p- 1 )Y Vl 



Sujlm, 



.r 



'L_i e --V- i ^(p;7)^ m - 



(2.19) 



For slow motion sources, we evaluate the asymptotic form of the Green's function up 
to O(Muj) as 



+0(MW)}e- i ^ t ~ r ~ t '\uryY l UmL(tt) + 0(r~ 2 ) 

'—— / dull + 2Mu[i(-Y.- + ln2M) + -s#n(w 



8=1 



r tr 2^(2/ + 

+i(ln \co\ + C)} + 0(MV)] e -^- r -*VO^^m(^) 
+0(r- 2 ), 



(2.20) 



where C is Euler's number. In deriving the above expression, we have used Eqs. fl2.16p 
and ( 2.17|) and the following expansion for q and 07 in Muj 

1 + ttMu + 0(M 2 u 2 ) 



Cl 



(21 + 1) 



1 1 

s=i s 



W) 



(2.21) 



Now we apply the formula pT| , |22 

u / dve luJV \nv + i / — e lu 

JO Jl V 
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— — sgn(uj) — i(\n \lu\ + C), 



(2.22) 



to Eq. (|2.20D , then we obtain 

1 ~ (-*y 



r f 2vr(2/ + 



— / d^{l - 2M(- S 1 + ln2M)| + 2M( £ d ve *- In, 



dt 2 



+2M( /°° — e^ v ) + 0(M 2 cu 2 )} 
xe— (*-»- *' V')'^)^') + 0(r- 2 ). 



(2.23) 



Here, we assume the no incoming radiation condition on the initial hypersurface so 
that we may take 
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lim e-^-^-t'hnv -> 0. 



(2.24) 



Thus we can make the following replacement 

1,00 dv ,, „, 



-IUJ 



i f 



dve iujv \nv= (^°° ffoe iaro In vj 



d_ 
~di 



(2.25) 



Inserting Eq.( |2.25| ) into Eq.( |2.23| ), we finally obtain the desired expression for the retarded 
Green's function 



i 



-i 



r tr 2 < 21 + 



I p ' 1 

-— / dw fl + 2M (X) - - In 2M 
1).. J 8=1 s 



eft 



+2M(^°° due*" In v) ^ + 0(M 



- part of 
r 



G (x, *') + 2M^- 2 / M ln (^l) + E -} G o m (^ - «, x, *' 

Zm s=l ^ 



+0(M 2 ) +0(r- 2 ), 



(2.26) 



where we defined the spherical harmonic expansion coefficient of the flat Green's function 



as 



G l m (x,x') = Y lm {n') / dn'Go(x,x')Y lm (n'). 



(2.27) 



As a result, we obtain the waveform generated by the effective source as 

72 „ „. 1-2 



h 



TT 



i 00 1 ~ d r v 1 

\P im j:^[nL-2{M pqL ^{t-r) + 2M— J du{ln(— )+g-} 



21 



xM pqL _ 2 (t - r - «)} - -p—^n aL _ 2 {tab( P S q ) b L-2(t - r) 

+0(M 2 )] +0(r- 2 ), (2.28) 

where -Pi^ is the transverse and traceless projection tensor and parentheses denote sym- 
metrization, and til is the tensor product of L radial unit vectors. Although this waveform 
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corresponds to (CI) by Blanchet | 23fl , M pqL _ 2 and S pq L_2 are the mass and current mul- 
tipole moments generated by the full nonlinear effective source f^ u . They take the same 
forms as in TJ], but here is used in place of t^ v in [I4"|, where 



;2.29) 



It is worthwhile to point out that In 2M in Eq. (|2.28|) can be removed by using the freedom 
to time translation. 



III. COMPARISON WITH THE PREVIOUS WORK 

A. waveform in the post-Minkowskian approximation 

By using the post-Minkowskian approximation, Blanchet obtained the general formula 
for gravitational waves including tail |[23|| , 

4 00 1 21 
h W = ~ Pi iPi 53 J\ [ n L-2U pqL - 2 (t - r) - j^n aL -2eab(pV q)bL - 2 (t - r)] + 0(r~ 2 ), (3.1) 

where the radiative mass moment and radiative current moment are defined as 

U [ l\u) = Mf\u) + 2GM ' dvMt +2 \u - v){\n[j) + «,} + 0(G 2 M 2 ), (3.2) 



vH\u) = S<£\u) + 2GM [°° dvS ( l +2 \u - v){ln(^) + k\) + 0(G 2 M 2 ). (3.3) 



P' 

Here the moments and Sl do not contain the nonlinear contribution outside the 
matter, P is a constant with temporal dimension, and H[ and k[ are defined as 



and 



_ if? 1 2Z 2 + 5Z + 4 
K '~£^ + Z(/ + l)(/ + 2) 



' _1 1 Z - 1 



(3.4) 



Z(Z + 1)' 



(3.5) 
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In black hole perturbation, the same tail corrections with multipole moments induced by 
linear perturbations have been obtained |[25|| . 

Equation ( [2. 28| ) does not apparently agree with Eqs.(|3.2|) and (|3.3|), because definitions 
of the moments {Ml, Sl} and {Ml, Sl} are different. We shall show below the equivalence 
between our expression and that of Blanchet [ 23| , by counting the contribution from 
nonlinear terms like M x Ml or M x Sl in the source f^ u . 



B. Contributions from the nonlinear sources 



In order to evaluate the tail of the waveform produced by the nonlinear sources in t^ v ' , 
it is enough to use the flat Green's function 

G (x,x') = J ^^(e-^/ i K^)^ m (fi)e^j i K / )F i ; i (^)^(^-^) 

-™j l (ur)Y lm (n)e^h l (ur')YC n (tt)e(T' " r) 



lm 



+e" 



(3.6) 



where ji and hi are the first kind of the spherical Bessel function and the spherical Hankel 
function, respectively. 

We consider the asymptotic form of the following retarded integral 



—i / duuje ^hiiur^FjfiL / r' 2 dr'ji(ujr' 

Jo 

for large r, 



Ik 



(3.7) 



where the hat denotes the symmetric and traceless part of tensor products and we defined 

F UJ as 



F^ = ±- fdte^F(t). 

Z7l J 



(3.8) 



Thus, the contribution from nonlinear sources can be evaluated by using the formula 



F(t - r) 



CO 

-2 l+1 limfV 



{l + n+ l)W(-k + I + 3 + 2n - A) 



n=0 



n\(2l + 2n + 2)\ 



) 



yV± {k F \ - r) + 0(r- 2 ) 
r 



(3.9) 



where we used 



^E 



[z/2) 



l+2n 



n=0 



n\T(l + n + 3/2)' 



(3.10) 



Since we consider both M x and M x Sl in the nonlinear source, we have to treat 
the following type of retarded integral 



□ " 



y \ nr 



E 



-i 



y 1 ^ + 2fTT 5i<a " nQ - 1> )^ F[t ~ r 



1 (q-3) 



(3.11) 



Using Eq. ( |3.9|) , we obtain 



^__ bi<ainQ ^ _ F(t _ r ) + (r- 2 ) 



(3.12) 



2(9+1) 

where, in the limit A — > of Eq. ( |3. 1 2|) , poles of gamma functions cancel out in total and 
the finite values are obtained. This is same with the expression (C2) obtained by Blanchet 



Applying Eqs.( |3.l2| ) to the nonlinear source ft*", we can evaluate its contribution to 
the waveform. Together with Eq.( 2.28|) , we obtain the waveform in the same form as 

Eq.Q. 



IV. DISCUSSION 



By considering the propagation of gravitational waves generated by slow motion 
sources in the Coulomb type potential, we have derived the formula Q2.28D for gravitational 
waves including tail. Our equation fl2.28|) has been shown to be same as Eqs. ([3.2|) - (|3.5|) 
in the previous work |23 . 



We would like to emphasize two points on Eq. (|2.28 ): First, in deriving Eq. (|2.28 ), 
spherical Coulomb functions are used, since we use the wave operator (|2.7| ) which take ac- 
count of the Coulomb-type potential M/r. As a consequence, hx(v/2M) appears naturally 
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in Eq. (|3.ip . This is in contrast with Blanchet and Damour's method, where an arbitrary 
constant with temporal dimension P appears in the form of \n(v/P). Our derivation 
shows that the main part of the tail, which needs the past history of the source only 
through lnt>, is produced by propagation in the Coulomb-type potential. The reason for 
saying the main part of the tail is as follows: Only the log term has a hereditary property 
expressed as the integral over the past history of the source, since the constants Ki and k[ 
represent merely instantaneous parts after performing the integral under the assumption 
that the source approaches static as the past infinity. Our method allows us to study the 
effect on the waveform formula by the modification of the Coulomb type potential near 
the material source. This effect might be important to make the comparison between 
theory and observation. We shall discuss this in future. 

The second point relates with physical application: At the starting point of our deriva- 
tion, the Fourier representation in frequency space has been used. Such a representation 
seems to simplify the calculation of gravitational waveforms from compact binaries in the 
quasi- circular orbit, since such a system can be described by a characteristic frequency. 
Applications to physical systems will be also done in the future. 
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